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The resonant-state expansion (RSE), a rigorous perturbative method developed in electrodynam¬ 
ics for non-dispersive optical systems is applied to media with an Ohm’s law dispersion, in which 
the frequency dependent part of the permittivity scales inversely with the frequency, corresponding 
to a frequency-independent conductivity. This dispersion has only a single pole at zero frequency, 
which is already present in the non-dispersive RSE, allowing to maintain not only the linearity of 
the eigenvalue problem of the RSE but also its size. Media which can be described by this dispersion 
over the relevant frequency range, such as optical glass or doped semiconductors, can be treated in 
the RSE without additional complexity. Results are presented using analytically solvable homoge¬ 
neous spheres, for doped silicon and BK7 glass, both for a perturbation of the system going from 
non-dispersive to dispersive media and the reverse, from dispersive to non-dispersive media. 

PACS numbers: 03.50.De, 42.25.-p, 03.65.Nk 


I. INTRODUCTION 

The electromagnetic spectrum of an open optical sys¬ 
tem is characterized by its resonances, having spectral 
positions and linewidths, corresponding to, respectively, 
the real (0) and imaginary (T) part of the complex eigen- 
frequencies u> = 0 — zT of the system. Finite linewidths 
(2T full width at half maximum) of resonances are typical 
for open systems due to energy leakage from the system 
to the outside, and are described by the quality factor 
Q = 0/2T. While the eigenmodes of resonators for a 
few highly symmetric geometries can be calculated ex¬ 
actly, simulating general system with high quality factors 
presents a significant challenge for the established simula¬ 
tion techniques such as the finite element method (FEM) 
or finite difference in time domain (FDTD) method. Re¬ 
cently, approximate approaches using resonant modes 
have been reported [ll-{5J. We have developed [g] a rigor¬ 
ous perturbation theory called resonant-state expansion 
(RSE) and subsequently applied it to one-dimensional 
(ID), 2D and 3D systems [7Hl0l|. demonstrating its abil¬ 
ity to accurately and efficiently calculate resonant states 
(RSs) - the eigenmodes of a perturbed open optical sys¬ 
tem, - using the spectrum of RSs of a simpler, unper¬ 
turbed one. 

Up to now, the RSE assumed a frequency independent 
permittivity. This is a significant limitation for the treat¬ 
ment of realistic systems which are all exhibiting disper¬ 
sion to some extend. In this paper we extend the RSE 
formulation from non-dispersive media to media exhibit¬ 
ing an Ohm’s law dispersion, i.e. having a permittivity 
consisting of a constant part plus a frequency dependent 
part scaling inversely with the frequency, as would be 
expected for a frequency-independent conductivity. This 
dispersion has only a single pole at zero frequency, which 
is already present in the non-dispersive RSE, allowing 
to maintain the linearity of the matrix eigenvalue prob¬ 
lem and its dimension. This can be understood as due 
to the property of this conductivity not to add memory 
to the system, and thus not to introduce additional ma¬ 


terial resonances a finite frequency which would act as 
an additional degree of freedom. Apart from conductive 
materials also other dispersive materials such as optical 
glass can be described by this permittivity over a rel¬ 
evant frequency range using an imaginary conductivity, 
and therefore can be treated in the RSE without addi¬ 
tional complexity. 

The paper is organized as follows. Sec. [II] outlines the 
RSE with an Ohm’s law dispersion. Examples illustrat¬ 
ing the method and comparing results with existing ana¬ 
lytic solutions are given in See. lIIIl Specifically, we show 
in Sec. lIII Al the doping of a silica sphere, and in Sec. IIII Bl 
an application to a dispersive glass sphere. 


II. RSE WITH OHM’S LAW DISPERSION 

Resonant states of an open optical system with a local 
time-independent relative permittivity (RP) tensor £fc(r) 
and relative permeability p = 1 are defined as the eigen- 
solutions of Maxwell’s wave equation, 

V x [V x E„(r)] = k 2 n e kn ( r)E„(r), (1) 

satisfying the outgoing wave boundary conditions. Here, 
k n = u n /c is the wave-vector eigenvalue of the RS num¬ 
bered by the index n, u> n is its complex eigenfrequency, 
and c is the speed of light in vacuum. In the present 
work, we consider the RP 

„ , , , *<r(r) 

£fc(r) = e(r) H - , (2) 

which includes additionally to the frequency-independent 
term e(r), a term with the conductivity <r(r), in accor¬ 
dance with Ohm’s law. 

A perturbed optical system is described by the wave 
equation 

V x [V x E(r)] = k 2 [e fc (r) + Ae fc (r)] E(r) (3) 
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with 

Ae fc (r) = Ae(r) + , (4) 

in which the RS wave vectors k and the electric fields 
E are modified as compared to k n and E n , respectively, 
due to a perturbation Aefc(r) inside the system. We treat 
this problem exactly by using the Green’s function (GF) 
Gfc(r,r') of the unperturbed Maxwell’s wave equation 
satisfying the outgoing wave boundary conditions, which 
yields a formal solution of Eq. (0: 


E(r) = —fc 2 / G fc (r,r')Ae fc (r')E(r')dr'. 


( 5 ) 


Substituting the spectral representation of the GF 

E„(r) ® E n (r') 


Gfe(r,r') = £■ 


2 k n (k — k n 


( 6 ) 


into the unperturbed Maxwell’s wave equation with a 
source term i<5(r — r') we obtain El the closure relation 


e(r) 


£ En(r) ® E„(r') = i<5(r - r') 


( 7 ) 


and the sum rule 


E„(r) ® E^r 7 ) _ 


( 8 ) 


both having the same form as in the non-dispersive 
case Q. Combining Eq. |6]) and Eq. (O yields Q 


Gfc(r,r') = Y 


E»(r) ® E n (r') 
2 k(k — k n ) 


( 9 ) 


This representation of the GF is then used in Eq. ([5]) 
leading to the following integral equation 


E (r) = £ 


E n (r) 

2 {k n - k) 


E„(r , )-fcAe fc (r , )E( r ')dr / . (10) 


Expanding the perturbed wave functions into the unper¬ 
turbed ones, 


E(i-) = ^c n E„(r), 


( 11 ) 


and using the linear independence of the basis functions 
E„(r), which holds for any finite basis size, we equate for 
each n the expansion coefficients at E„(r) in both sides of 
Eq.GDJ. For the perturbation in the form of Eq. (j4j , this 
results in a linear generalized matrix eigenvalue problem 

^ ) (k V n m + iS nrn )c rn = 2(/c n k)c n , (12) 

m 

in which 

V nm = j E„(r) • Ae(r)E m (r) dr , (13) 

S nm = j E„(r) • A<r(r)E m (r)dr . (14) 


Without perturbation in the conductivity, S nm vanish, 
and Eq. m reduces to the previously known matrix 
equation of the RSE, which can be reduced to a sym¬ 
metric matrix eigenvalue problem [see Eqs. (14) and (15) 
in Ref.0- All modes appearing in Eqs. © (fl4l) are nor¬ 
malized 0] according to 


1 = 


E„(r) 


d(k 2 e k ( r)) 


>v 


d(k 2 ) 


E„(r)dr 


k—k n 


J Sv 


Er, 


d 9E„ 


<9E„ dE r 


ds dr 


dr ds 


(15) 

dS 


where the first integral is taken over an arbitrary simply 
connected volume V enclosing the inhomogeneity of the 
system and the center of the spherical coordinates used, 
and the second integral is taken over its surface Sy with 
d/ds being the gradient normal to the surface. This nor¬ 
malization is required M for the validity of the spectral 
representation Eq. ([9]). 


III. APPLICATION TO SPHERICAL 3D 
SYSTEMS WITH SCALAR PERMITTIVITY 


In this section we show applications of the simple dis¬ 
persive RSE to 3D systems described by a scalar RP, i.e. 
£fe(r) = 1 £fc(r), and Aefc(r) = lAefe(r). As unperturbed 
system we use a homogeneous sphere of radius R and 
efc = e + ia/k in vacuum, so that 


( e k for r < R , 
( 1 elsewhere. 


(16) 


We treat transverse magnetic (TM) and transverse 
electric (TE) modes, while longitudinal electric (LE) 
modes which all have k = 0 are not present in the dis¬ 
persive basis due to the divergence of e k at k = 0 and 
vanishing normalization coefficients. 

Owing to the spherical symmetry, the modes can be 
classified according to the orbital quantum number l and 
magnetic quantum number to. The modes E„(r) of this 
system are of the same analytical form as in Ref. ly. The 
resonant wave numbers k n are given by the solutions of 
known EJ secular equations, see for example Eqs. (30) 
and (31) in Ref. For a dispersive e k the normaliza¬ 
tion coefficients A TE and A™ given in Ref. S3 are mod¬ 
ified El 1° 


A 


TE 

i 


1(1 + 1)-R 3 (efc — 1 + r)Bi) 


(17) 


a tm 


2 e k 


l(l + l)R 3 (e k -l) 


/ 


( Ji-i(z) 

V M z ) 



+ Cfc 


1(1 + 1 ) 


■vCi 


(18) 
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where z = y/ckkR , and the additional terms in 
malization coefficients due to the dispersion are 

the nor- 
given by 

(a) 


1 d{k 2 e k ) 

V e k d(k 2 ) 

(19) 

Ph 

Pd 


Bi= e k \i , 2 

V 3i(z) J 

(20) 

(b) 

Ci = 

/ 0 i + 1 3?+ i( 2r ) jl+ 2 (z)\ 

(21) 

£ 

¥ 

i 

e fc — 1 \ * 2 jf(z) jl(z) )' 


20 


where all the quantities above are evaluated at k = k n . 


A. Doping of a silicon sphere 

Here we show as example a homogeneous n-doping of 
a silicon sphere at room temperature. The unperturbed 
system is a sphere of undoped silicon with e = esi = 
(3.42) 2 . The intrinsic conductivity of about 10” 4 S/m is 
neglected and we use a = 0. Introducing an n-doping 
with a density of 10 14 /cm 3 leads to a conductivity of 
crgi = 2.3S/m [l3|. The Drude scattering time is about 
200 fs for this doping density [l3| , dominated by phonon 
scattering, such that the Ohm’s law description is valid 
within 10% error for frequencies below 200 GHz, corre¬ 
sponding to \kR\ < 21 for R = 5 mm. The resulting 
perturbation is Act = crsi/(eoc) = 4.334/1? with the vac¬ 
uum permittivity eo, while Ae = 0. The resulting RP of 
the doped silicon is shown in Fig.QJa) as function of kR. 

For perturbations retaining the spherically symmetry, 
RSs of different angular quantum numbers l , to and dif¬ 
ferent transverse polarizations are not mixed, and are 
degenerate in to. The procedure used for the RSE cal¬ 
culation follows that described in Ref. ES Namely, we 
choose the basis of RSs for the RSE in such a way that 
for given Z, m and TE or TM polarization, we select all 
RSs with \k n \ < Zc max (iV) using a maximum wave vector 
kmax(N) chosen to select N RSs. 

The resulting RS wave vectors are given in Fig. [lib) 
for TE and Fig. [lie) for TM states, for the unperturbed 
system (Zc n ) and for the perturbed system (Zc). The per¬ 
turbed wave vectors calculated using the RSE with N = 
200 are compared with the exact result fc( exact ) from the 
secular equations. The relative error S = |fc/fc( exact ) — l| 
is shown in Fig.[lld) for both TE and TM polarizations. 
We find that as we increase N, S decreases proportional 
to N “T similar to the findings for the non-dispersive RSE 
an d values below 10 -8 are reached for N = 200. 

Interestingly, the Ohms law dispersion gives rise to 
nearly constant imaginary part of the whispering gallery 
modes (WGMs) for large kR, which is due to the con¬ 
tribution of the absorption in the medium. In detail, 
k = k 0 /y/ek = ko/y/e+lajk « k 0 /\/e(l - icr/(2ke), for 
|cr/(Zce)| <C 1, where fco is the peak wave vector of the 
dissipated wavepacket, resonant to the frequency of the 
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FIG. 1: RSE results for the doping of a silicon sphere, as 
function of the real part of the dimensionless wave number kR. 
(a) Real and imaginary part of the RP t k of doped silicon. 
(b,c) RS wave numbers for TE (b) and TM (c) polarization 
and N = 200. Shown are exact unperturbed (circles), exact 
perturbed (squares), and RSE-perturbed (crosses) data. The 
symbol color indicates l as labeled, (d) Relative error 5 of the 
RSE-perturbed RS wave numbers k for N = 50,100, and 200, 
for TE and TM modes with Im (kR) < 1, as labeled. 


given RS: k ~ ko/y/e. Then the imaginary part of the 
RS wavenumber is approximately —<r/(2e), which has a 
value of about 0.18/i? for the doped silicon considered. 
This value is close to the imaginary part of WGMs seen 
in Figs. [l]b) and (c) for kR > 5. The Fabry-Perot modes 
(FPMs) have additional losses due to the finite reflection 
at the silicon to vacuum interface, resulting in an ad- 
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FIG. 2: RSE results for the un-doping of a silicon sphere, 
as function of the real part of the dimensionless wave number 
kR. (a) RS wave numbers for both TE (filled symbols) and 
TM (open symbols) polarization and N = 200. Shown are 
exact unperturbed (circles), exact perturbed (squares), and 
RSE-perturbed (crosses) data. The symbol color indicates l 
as in Fig. |T| (b) Relative error 5 of the RSE-perturbed RS 
wave numbers k for N = 50,100, and 200, for TE and TM 
modes with Im (kR) < 1, as labeled. 


ditional imaginary part of about 0.09/i?. Leaky inodes 
(LM) are mostly outside the sphere and are therefore less 
influenced by the absorption. 

We now exemplify that the RSE can also be used with 
a basis from a dispersive system, which is perturbed to 
a non-dispersive one. We use the //-doped silicon sphere 
as unperturbed system and perturb it by removing the 
doping, i.e. A a = —4.334/1?. The resulting RS wave 
numbers and their relative error are given in Fig. [2] for 
different basis sizes N. We find that the wave numbers 
are well reproduced, and that the relative errors are scal¬ 
ing equivalently to the case of a non-dispersive basis. 


B. Dispersive glass sphere 

The form of the dispersion, representing Ohm’s law 
with a frequency independent conductivity, is specifically 
suited to describe conductors at frequencies below the 



FIG. 3: RP n 2 of Schott BK7 glass (solid line), and fitted 
Ohm’s law type dispersion (dashed line). 

charge carrier scattering rate. However, the mathemat¬ 
ical treatment is in general allowing for complex e and 
<t, which can be used to model other types of materi¬ 
als. As example we use a technical glass in the optical 
frequency range, specifically the borosilicate glass Schott 
BK7, which has a refractive index n r in the optical fre¬ 
quency range described by the Sellmeier expression 0 

2 „ 1.03961212A 2 

n r — 1 + T 2 -conn R0SC7-2 (22) 

A z — 6000.69867 nnr 

0.231792344A 2 | 1.01046945A 2 

+ A 2 - 20017.9144 nm 2 + A 2 - 103.560653//m 2 ' 

The resulting RP is real and is shown in Fig. [3] It is fitted 
by efc using the parameters e and a over the wavelength 
range from 1 //m to 1.8 //in. The resulting e kl which is lin¬ 
ear in wavelength A = 2n/k 7 is also shown in Fig.[3j and is 
given by e = e g = 2.30926 and a = cr g = 0.232414i///m. 
The deviation of the fit is below 0.001 over the fitted 
range. One has to keep in mind that a complex conduc¬ 
tivity results in an unphysical dispersion, which does not 
have the property e k = e*_ k , following from the causality 
principle [15l | . The latter requires in particular that the 
conductivity has to be real. However, this does not af¬ 
fect the suitability of the RSE to describe the RSs close 
to the fitted k range, and we have evaluated the error 
of the fit on a selected segment of k. We note that the 
value of the permittivity at the generally complex k n will 
deviate to some extend from a physical RP which obeys 
Kramers-Kronig relationships. 

We first show RSE results introducing dispersion as 
perturbation of a non-dispersive system. The unper¬ 
turbed system is given by a sphere with e = e g , cr = 0 
and R, = 7 //m. We consider a homogeneous perturba¬ 
tion over the sphere, having Act = cr g , Ae = 0, so that 
the perturbed system has the fitted BK7 dispersion. We 
choose the basis of RSs for the RSE in the same way as 
in Sec. lIII Al The results are given in Fig. [I] In Fig.|U(a) 
the refractive index of the unperturbed and perturbed 
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FIG. 4: RSE results for a homogeneous sphere with a non- 
dispersive unperturbed and a dispersive perturbed medium, 
as function of the real part of the dimensionles wave number 
kR. (a) Refractive index of the unperturbed (dotted line) 
and perturbed (dashed line) medium, together with the BK7 
Sellmeier dispersion (solid line) given by Eq. (1221) . (b,c) RS 
wave numbers for (TE,TM) polarization and N = 200. Shown 
are exact unperturbed (open circles), exact perturbed (open 
squares), and RSE-perturbed (crosses) data, (d) Relative er¬ 
ror S of the RSE-perturbed RS wave numbers k for N = 
50,100, and 200, for TE and TM modes with Im(fcit) < 0.1, 
as labeled. For the N = 200 data the relative error to the 
exact RS wave numbers for the Sellmeier dispersion is also 
shown. 



Re(kR) 


FIG. 5: RSE results for a homogeneous sphere, with a dis¬ 
persive unperturbed and a non-dispersive perturbed medium, 
as function of the real part of the dimensionless wave number 
kR. (a) RS wave numbers for both TE and TM polarization 
and N = 200. Shown are exact unperturbed (open circles), 
exact perturbed (open squares), and RSE-perturbed (crosses) 
data, (b) Relative error $ of the RSE-perturbed RS wave 
numbers k for N = 50, 100, and 200, for TE and TM modes 
with lm(kR) <0.1 , as labeled. 


system is shown as function of kR, and compared with 
the BK7 Sellmeier data. We can see that for the chosen 
sphere radius, the fitted wavelength range corresponds 
to 22 < kR < 55. The RS wave vectors are given in 
Fig-ffib) for TE and Fig.@Jc) for TM states, for the un¬ 
perturbed system ( k n ) and for the perturbed system, and 
the relative error S is shown in Fig.QJd). We again find 
that as we increase N, S decreases proportional to N~ 3 , 
reaching values in the 10~ 8 range for N = 200. This is 
actually smaller than the relative error due to the Ohm’s 
law approximation of the Sellmeier dispersion, which is 
in the 10 -4 — 10 -5 range, as also shown in Fig.QJd). For 
practical applications, the simple dispersive RSE is thus 
limited by the Ohm’s law approximation of the dispersion 
over a given spectral range, but the achievable accuracy 
in the 10 -4 range can be sufficient for many applications. 

Also for this example we now verify that one can also 
use a dispersive basis to treat a dispersive perturbation. 
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We revert the dispersive perturbation, so that the unper¬ 
turbed system is given by a sphere with e = e g , a = a g 
and the perturbation is Act = — cr g , Ae = 0. The result¬ 
ing RS wave numbers and their relative error are given in 
Fig. [5] for different basis sizes N. We find that the wave 
numbers are well reproduced, and that the relative errors 
are scaling equivalently to the case of a non-dispersive 
basis shown in Fig. [I] 

IV. SUMMARY 

In this work we have extended the RSE to media hav¬ 
ing a simple dispersion equivalent to Ohm’s law, i.e. a 
frequency independent conductivity. This dispersion has 
a single pole at zero frequency and does not introduce 
additional dynamic degrees of freedom. This property al¬ 
lows to keep the simplicity of the RSE formulation, with 
the only change to move from a normal matrix eigen¬ 
value problem to a generalized matrix eigenvalue prob¬ 
lem, keeping the matrix size the same. With suited al¬ 
gorithms, this results typically only in a doubling of the 
computational time, therefore retaining the advantage of 


the RSE in computational efficiency discussed in Ref. 0. 
The Ohm’s law is a good approximation for conductive 
materials at frequencies well below the carrier scattering 
rate, which is typically in the terahertz frequency range 
at room temperature, and we show as example doped 
silicon. Furthermore, we have shown that by choosing 
an imaginary conductivity, this dispersion can also de¬ 
scribe the dispersion of the real part of the permittivity 
of optical glass in the normal dispersion range, and the 
RSE is working equally well despite the dispersion not 
respecting causality. 
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